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Abstract. We use molecular dynamic to simulate the directional growth of binary 
mixtures, our results compare very well with analitical and experimental results. 
This opens up the possibility to probe growth situations which are difficult to reach 
experimentally, being an important tool for further experimental and theoretical 
developments in the area of crystal growth. 

1 Introduction 

The rapid expansion of the use of high quahty crystahne materials in optical 
and eletronic devices has strongly stimulated research, both theoretical and 
experimental, on dynamics of crystalization. Computer simulation has played 
an important role on the development and understanding of crystal growth. 
During growth the solid-fluid interface can display several interesting phe- 
nomena like segregation, dynamical instabilities and pattern formation | ^ . 

A crystal can grow from the adjacent fluid by different mechanisms, de- 
pending on the structure of the interface (rough or smooth), material purity, 
growth rates, temperature gradients and so on. For a crystal to grow atoms 
or molecules must be transported from the fluid towards the interface with a 
non-zero sticking probability, besides that, the latent heat generated at the 
interface as well as the solute excess segregated must be carried away. These 
requirements can be met in a controlled way by puting the sample in an appro- 
priated furnace submmited to a temperature gradient and pulling it with con- 
stant speed toward the colder region, results of such experiments have been 
compared with results of two dimensional models. Nijmeijr and Landauj^ re- 
ported molecular dynamic (MD) simulation on laser heated pedestal growth 
of fibers. Previous simulations consisted of kinetic models using Monte Carlo 
techniques and numerical solution of transport equations B. 



2 Simulation 

In this work we use MD to simulate the solidification of a two component 
system of a solvent (particles of type a) and solute (particles of type b) in- 
terating via a Lennard-Jones (LJ) potential. By tunning the LJ parameters 
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we set the struture of the interface (rough or smooth) and the segregation 
coefficient. We also simulate morphological instabilities, the planar inteface 
becoming cellular and eventually dendritic (MuUins-sekerka instability.). The 
simulation is carried out with all particles interacting through the LJ poten- 
tial 



'/'^J(''^J) 
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The indexes i and j stands for particles in the positions and rj re- 
spectively and ri_j = \ri — There are three types of interactions, solvent- 
solvent, solvent-solute and solute-solute. In each case the LJ parameters are 
labeled as (a, a), (a, 6) and {h,b). distance r, time t and temperature T are 
measured in units of Oa.a, <ya,a('ma/ ^a.a)^ and kB/^a.a- initially we distribute 
N = Tlx X riz particles over the two dimensional volume x L^. We as- 
sume periodic boundary conditions in the x direction. In the z direction we 
divide the system in to two distinct regions, a solid and a fluid one. In the 
solid region particles stand initially in their equilibrium position in a total of 
nx X particles, on the fluid region we distribute the rest of the particles 
with the density chosen to be p = O.bcr^^, randomly distributed in a triangu- 
lar lattice and slightly dislocated from their equilibrium position. We impose 
a temperature gradient along the z direction using a velocity renormalization 
approach. The system is divided in two regions, one cold which is set to T = 
and the other to T = higher than the melting temperature which is 
Tm = 0.4 in our units, we let the system evolve in time for time steps of 
size 5t until it reachs equilibrium. Once the equilibrium is reached we start 
pulling the system in the -\-z direction, at a pulling velocity Vp. The —z^ax 
layer works as a particle source, maintaining a constant flow to the material. 
The basic experimental set up is shown in Fig. |l| 

In Fig. H we show an example of solute segregation during directional 
growth o fthe binary mixture: caprolactane as solvent and methyl-blue as 
solute. In the top part of Fig. || is shown an image of the crystal (left side) 
and melt (right side), with maximum concentration of methyl-blue at the 
melt side of the interface. From the gray level of the image we obtain the 
methyl-blue concentration profile across the sample, which in the melt, decays 
exponetially as a function of the distance from the interface (bottom part of 
Fig. H). Fig. m shows the simulation results done in a system initially with 
27 X 30 particles, solute concentration of 5%, pulling velocity Vp = 0.004, and 
LJ parameters given by : ta^b = 0.5, €b,b = 0.1, (Ja^b = 1-0, (Jbfi — 1-0 and 
mb = 1.0. 

After averaging over many runs to improve statistics one obtains the 
steady-state solute concentration profile represented as data points with error 
bars. 



Molecular Dynamic Simulation of Directional Crystal Growth 



3 




Fig. 1. Basic experimental set up for directional growth. The interface motion is 
visualized using an optical microscope. Tm is the melting temperature of the mixture 
and Vp is the pulling velocity. 
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Fig. 2. Directional solidification of a binary mixture showing the solute segregation 
at the interface. The experimental result is for the mixture caprolactane (solvent) 
and methyl-blue (solute). 
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Fig. 3. Our simulation results using the Lennard- Jones potential describes in detail 
the experiment. 



It is shown in Fig. ^ a theoretical resuh obtained from the diffusion equa- 
tion for the solute concentration in the fluid phase jij 

where D is the solute diffusion coefficient in the fluid phase and Vf is the 
velocity of the system of reference. The diffusion equation is supplemented 
with the boundary conditions: 



c/ = Co at z ^ GO , (3a) 

{l^ K)VfCf = -D-^ at z = 0. (3b) 

The second condition is just the mass conservation at the interface and 
K is the ratio between the solute and solvent concentrations. 



3 Cellular Instability 



With the LJ parameters chosen above the morphological instabilities were 
inhibited during the cellular growth. We can stimulate instabilities by tunning 
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the LJ parameters, particularly by slightly varying the equilibrium position 
between particles a and b, i.e., changing aa.b- Fig. ^show the effect of varying 
this parameter. 
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Fig. 4. The figure shows the concentration profile found in our simulation. It is to 
be compared to Fig. |^ 



For CTa,h — 1 the interface is rough, becoming cellular and dendritic for 
Ca,6 > 1. As shouls be expected a lower cTa.t sets the interface smoother than 
larger values indicating that the system is very sensitive to geometric factores. 
(We have also simulated some experiments fixing the value of aa,b and varying 
Ca.b and eb^b, the results did not show any qualitative change in their behavior. 
Raising (7a, 6 breaks the hexagonal structure o fthe solid. The energy involved 
in such distortion is so large that the interface segragate solute particles, 
bonds (a, a) and (6,6) are preferible than (a, 6). Once particles of the type 
6 concentrate at the interface its melting point is lowered and channels of 6 
particles are formed inside the crystal. Fig. || Shows this effect. 

We are now studying the stability 1^ of those interfaces. We start with a 
smoth interface then, we introduce a perturbation at the interface. It has been 
shown that near the bifurcation from planar to cellular the time evolution 
of the most instable Fourier mode of the perturbation can be described by 
a third order Landau-amplitude equation. This however is the subject of a 
new research and will be soon published elsewhere. 
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